Electromagnetic Casimir Energies of Semi-Infinite Planes 
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Using recently developed techniques based on scattering theory, we find the electromagnetic 
Casimir energy for geometries involving semi-infinite planes, a case that is of particular interest 
in the design of microelectromechanical devices. We obtain both approximate analytic formulae 
and exact results requiring only modest numerical computation. Using these results, we analyze the 
effects of edges and orientation on the Casimir energy. We also demonstrate the accuracy, simplicity, 
and utility of our approximation scheme, which is based on a multiple reflection expansion. 

PACS numbers: 

Introduction and Method The Casimir energy has been most commonly studied using techniques based on the 
original calculation for infinite parallel planes Recently developed techniques @, 0] have made it possible to 
■ calculate the electromagnetic Casimir force between objects of arbitrary shape and electromagnetic response. In recent 
, applications of these methods, the interaction energy of a semi-infinite plane with an infinite plane was analyzed in 
two different ways: First, with the half-plane considered as the limit of a parabolic cylinder of zero curvature in 
which case the exact energy can be computed numerically. Second, with the semi-infinite plane taken as the limit of 
a wedge of zero opening angle in which case it is convenient to consider a multiple reflection expansion for the 
I , energy. Remarkably, the analytic formulae obtained by keeping the lowest few orders in the reflection expansion give 
very good agreement with the full numerical result. (This approach can also be extended to interactions of more than 
two bodies [6j.) We study the Casimir interaction of semi-infinite planes by contrasting these two methods. This 
^ ■ problem is applicable to the design of microelectromechanical devices Q and has been of recent theoretical interest 
as well d, Q . We present analytic formulae obtained through the multiple-reflection approximation and use the exact 
numerical calculation to obtain a concrete measure of the accuracy of the approximations involved. This system 
-H , provides an ideal environment in which to study the effects of edges and orientation. 

In the scattering theory approach [3, [3| , the Casimir interaction energy is expressed in terms of the scattering T- 
matrices, also known as scattering amplitudes, for each object individually. These matrices incorporate the material 
characteristics of each object individually, while the objects' relative positions and orientations are described by 
universal translation matrices. These matrices connect the bases of wavefunctions, centered on each object, in which 
1 the scattering amplitudes are calculated. For the case of two objects that are translationally invariant in the z- 
\ direction, we consider the energy per unit length 
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where the Tj give the T-matrices for each object, while the Uij = L{j i give the translation matrices from one object 
to another. Here k = in = oj/c is the magnitude of the wave vector for each possible fluctuation and k z is its 
^-component. We can interpret this formula in terms of the propagation of electromagnetic fluctuations between 
the objects: The T-matriccs describe the reflection of fluctuating fields from a single object, while the ^/-matrices 
propagate these fluctuations from one object to another. The determinant then combines all possible reflections 
among the objects. By evaluating the integral and determinant numerically, we can use Eq. fTJ) to find exact results 
for the Casimir interaction energy. However, this expression also allows for a systematic approximation, in which we 
can obtain simple analytic results. Letting Af = T^UyiT-JUix-, we can convert the log-determinant to a trace-log, which 
we then expand as a Taylor series, 
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In this expansion, the n th term then gives the contribution from n reflections (back and forth) between the two objects. 
These contributions typically fall at least as fast as n~( D+1 \ where D is the dimension of space [Bl [Toj. 

We can realize the half-plane geometry as the limit of either a wedge [B| with zero opening angle or a parabolic 
cylinder Q with zero radius of curvature. In the former case the T-matrix is expressed in a basis indexed by 
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continuous imaginary angular momentum A, while in the latter case the T-matrix is given in a basis with discrete 
channels v = 0,1,2.... Although we can compute either the full determinant or the multiple reflection expansion 
using either basis, the wedge basis is better suited to the reflection expansion, because the associated translation 
matrix elements are simpler and easier to handle analytically, while the parabolic cylinder basis is better suited to 
the full determinant calculation, because it is easier to calculate the determinant when the matrix involves discrete 
rather than continuous indices. 

Two Half-Planes We consider two half-planes, and restrict our attention to configurations that are translation 
invariant in the z direction. We introduce a translation d v in the y direction and d x in the x direction, and allow 
the upper and lower half-planes to rotate around their edges away from the y-axis by angles 9 and 9 respectively, 
as shown in the left panel of Fig. [TJ This description is redundant — different parameter choices that lead to the 
same physical configuration, a property we use to check our calculations. The numerical convergence of physically 
equivalent configurations is not necessarily equivalent, however. For example, in the case of 9 = 8 = 0, when both 
d x and d y increase, the Casimir interaction energy decreases, since the planes are becoming further apart. For the 
scattering bases we choose, however, in dy this effect appears directly through a decaying exponential, while in d x it 
appears through the cancellation of an oscillating integrand. As a result, we need to maintain d y > 0, while we can 
consider either sign of d x . 




FIG. 1: Geometry for two half-planes (left panel) and a half-plane opposite a plane (right panel). 
For a wedge with zero opening angle, the T-matrix becomes [H| 

T D / N - zn - (3) 
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where the channels are labeled by Dirichlet and Neumann boundary conditions, corresponding to the two polarizations 
of the electromagnetic waves, and by the continuous index A, the imaginary angular momentum. We can then express 
the matrix Af in the wedge basis as 

poo 

Af D/N xy(d x ,d y ,8,9)= / d\''T^ N U^(d x ,d v ,e^)T^U^(-d x ,d y ,-e,-e), (4) 



where n and k z have been suppressed because all the matrices involved are diagonal in these parameters. Without 
loss of generality, d x can be set zero, in which case the translation matrix IA is given by 

r°° j 

U{ D x ,(O,d v ,6,0)=2 dk x — cosh(\{<j,- 9)) cosh(\' (0 + 9))* e lk v d y , (5) 
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for the polarization corresponding to Dirichlet boundary conditions, where k y = ^/k 2 — k 2 — k 2 = i^J k 2 + k 2 , + fcS? 
and (f> ~ sin" 1 (k x /i^Jn 2 + k^j. For Neumann boundary conditions, the hyperbolic cosine functions in Eq. ((5]) are 
replaced by hyperbolic sine functions. 
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In the parabolic cylinder basis, Af becomes 

00 

Af vl/ '(d x ,d y ,6,9) = ^ Tx u U vv >t(d x ,dy,6,d)T-z v >:Uv» v >{-d x ,dy,-6,-6) , (6) 
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where T av = —\j^v\ is the T-matrix for the half-plane a = 1, 2 in the scattering channel v = 0, 1, 2, 3 . . ., with even 
v corresponding to polarizations obeying Dirichlet boundary conditions and odd v corresponding to polarizations 
obeying Neumann boundary conditions, and the translation matrix is 

(tan^fftan^V' 

U^(d x ,dy,6,9)= dk x -± JLA L e ik*d* e ik v d (?) 
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where <f> is defined as before. 

We first evaluate the Casimir energy analytically to first order in the reflection expansion. Using Eqs. ([3]), ((H), and 
([5]). the electromagnetic Casimir interaction energy becomes 

— = \r-fr ( -+4:Csc9csc9 + 4:(6csc 2 6-9csc 2 9)csc(9-9) ) +••• , (8) 
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where without loss of generality we have taken d x — and the dots represent higher reflections. Since the contribution 
at each reflection order is independent of the scattering basis in which it is computed, this result can be obtained using 
cither the wedge or parabolic cylinder basis; however, the wedge basis is more convenient because it yields integrals 
rather than sums. In the parabolic cylinder basis, one can obtain the same results by writing the sum over v as a 
geometric series in tan . 

We are particularly interested in the case of parallel, overlapping half-planes. In this case, it is more convenient to 
fix 6 = 9 = tt/2 and parameterize the configuration by d x and d y . Then positive d x describes the width of the overlap 
region, while negative d x describes a horizontal displacement of the edges away from each other. The exact Casimir 
interaction energy and the approximation to one reflection for this case are shown in Fig. [2] At large d x , the graph of 

^£ as a function of d x /d y becomes asymptotic to a straight line with slope — §^tti consistent with the standard result 
for parallel planes: The energy is linear in the exposed area and inversely proportional to the cube of the separation 
distance, 

hcL 720c23 8n 2 d y ' [ ' 

This result forms the basis for the proximity force approximation (PFA), which sums over infinitesimal segments 



treated as locally parallel planes [Hj . The y- intercept of the asymptote gives the correction — beyond PFA — due to 
the interaction between each of the two edges and an infinite plane. It represents a constant shift of the parallel plane 
result in the limit of large d Xl where the edges are far from each other. This correction is positive — it makes the overall 
Casimir interaction energy less negative — reflecting the suppression of quantum fluctuations in the neighborhood of 
the sharp edge. The dashed line in Fig. [5] shows the result obtained from the first reflection, 
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which is Eq. ((HJ specialized to the geometry of overlapping parallel planes and expressed in terms of d x and d y . Here 
the logarithm gives the arctangent in the appropriate quadrant. For large d x , this result approaches the straight line 

— = (11) 
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This equation represents the first reflection approximation to the full result for parallel planes, which is obtained by 
replacing £(4) in Eq. (|9]) with the first term in the expansion of the zeta function, [To| 
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We can see from Fig.[2]that the first reflection already gives an excellent approximation to the energy of the overlapping 
planes. 

Note that two half-planes also exert a lateral force on each other (see Ref. [l2| for a gedanken experiment based on 
this point). For positive d x (when the planes overlap), and to leading order, this force can be obtained from the PFA, 
Eq. or, to one reflection, from Eq. (jlip . The y- intercept, which quantifies the interaction of each edge with an 
infinite plane, does not contribute to the lateral force, since it is independent of the distance between the two edges, 
d x . However, the interaction between the edges makes a contribution. The analytic formula of Eq. (fTU|) gives 



^=8^-6^( 1 + °(4M))+- d x >0. (13) 

The first term is simply PFA in the first reflection, while the second term gives a negative correction, representing a 
suppression of the force due to the interaction between the two edges. For negative d x , Eq. (|10l) gives 

We see that in the first reflection, the leading term of the lateral force for d x < is the same as the sublcading term 
in the force for d x > 0. One can use an argument based on Babinet's principle to show that this result holds for 



arbitrary d x /d y . Applications of Babinet's principle to the Casimir energy are studied extensively in Ref. [13 1. 

Half-Plane and Infinite Plane The edge correction can also be obtained by considering the limit in which a half- 
plane becomes parallel to an infinite plane. We consider a half-plane separated by a distance d from and tilted at an 
angle 9 from the perpendicular to an infinite plane, as shown in the right panel of Fig. [TJ As 9 — > ir/2, the Casimir 
energy per unit length diverges, since the energy is becoming proportional to the area, but we can parameterize this 
divergence as 

HcL cosfld 2 ' V 1 

where d gives the separation between the edge of the half-plane and the infinite plane, and 9 gives the angle between 
the half-plane and the axis normal to the full plane. Then for 9 — s- n/2, where the planes become parallel, we have 
c(9 —> 7r/2) = C||/2 + (9 — 7r/2) c odgo , where cm = 7r 2 /720 is the standard result for parallel planes. Then the slope 
Cod g o gives the correction due to the interaction between the edge of the half-plane and the infinite plane, in the limit 
where the half-plane is parallel to the infinite plane. Doubling this result gives the edge correction for the case of 
overlapping planes, since there we have two edge-plane interactions. 

The calculation for a half-plane opposite an infinite plane proceeds analogously to the case of two planes, and has 
been done in detail in Rcfs. [4j,[5[- In Fig. [3] we show the coefficient c(9) using both the exact numerical calculation 
in the parabolic cylinder basis jj| and the approximation to two reflections in the wedge basis [H[ , which is 

secfl- 1 \\ + csc 3 9 sec 9 (29 -sin 29) ) +••• , (16) 
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where the dots indicate corrections from higher reflections (three reflections or more). 

The comparison of the exact numerical calculation in the parabolic cylinder basis to the reflection expansion in 
the wedge basis shows both the strengths and weaknesses of this approximation. The expansion to two reflections 
for the energy of a tilted half-plane opposite a plane captures the orientation dependence c(9) well. The variation 
in this quantity is small because of the near-cancellation of the Dirichlet and Neumann contributions. However, the 
relative error in the edge correction, that is, the error in the derivative of c(9) as 9 — > 7r/2, is large because while 
the neglected contribution from the higher reflections is small, it varies rapidly as the planes become parallel. For 
9 = it/2, the reflection expansion to order n gives the parallel-plane result with £(4) truncated to the first n terms 
in its series representation, but this error decreases as the angle deviates from tt/2. This behavior can be anticipated 
from a geometric optics point of view, since as the planes become nonparallcl, higher-order reflections must propagate 
further as they reflect between the planes. The analogous problem arises in extracting the edge correction from the 
overlapping planes calculation, since again the asymptotic result as d x — > oo captures only the leading term in the 
zeta-function series. 

Thermal Corrections We can also apply these results to find the free energy in a system at temperature T. Then 
the integral over k is replaced by a sum over Matsubara frequencies k„ = 2imT '/he, where n = 0, 1, 2, . . . and the n = 
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FIG. 2: Electromagnetic Casimir interaction energy per unit length for overlapping planes as a function of horizontal displace- 
ment, in units of the vertical separation d y . Circles are obtained from the exact calculation in parabolic cylinder coordinates. 
The solid line connecting them is a rational function fit to guide the eye. The dashed line gives the analytic formula obtained 
by considering only the first reflection, while the solid straight line gives the standard parallel plane result for the overlap area, 
plus edge corrections. 
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FIG. 3: The orientation dependence c(8) for the electromagnetic Casimir energy of a half-plane opposite a plane, as a function 
of the tilt angle 6. Circles are obtained from the exact calculation in parabolic cylinder coordinates, while the solid line is 
obtained from the analytic expansion through two reflections, Eq. (|16[) . 



contribution is counted with a weight of 1/2. Again, we keep only the first reflection. For the case of a half-plane 
opposite a plane, we obtain the free energy for each polarization 
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where the dots represent higher reflections. The first term in parentheses cancels in the sum over polarizations, giving 
for electromagnetism 



EM 



hcL 



1 1 1 , d 11 

. cotli = 

32tt 2 cos X T d A T 32n 2 d 2 cos 9 



1 



3 \Xt 



1 



45 V At 



r/ (1 



(18) 



where At = he/ {iirksT) is the thermal wavelength. The leading correction at small temperature is proportional to 
T 2 , but this term is independent of distance and thus does not contribute to the force, for which the leading correction 
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starts at order T 4 . We note that in the case of a scalar field with either Dirichlet or Neumann boundary conditions, 
however, the first term in parentheses in Eq. (|17p yields the contribution 
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where E(x) is the elliptic function. This expression shows nonanalytic behavior as T — > proportional to T 3 lnT, 
similar to what was observed in Ref. Q. 
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